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^ ■ Abstract 

I \ Particles are either bosons or fermions, which obey Bose-Einstein (BE) 

Q \ or Fermi-Dirac statistics, respectively. In the hadronization process the 

1-^ ' effects of the former kind are the most important ones, since mesons are the 

^ ■ most common products in high-energy collisions. In the phenomenological 

^ . Lund Model the BE effect is approximated by a semiclassical momentum 

dependent correlation function, which effectively acts as an attractive force 
between two mesons. However, the Lund Model does also provide a space- 
time picture of the hadronization process. Based on this, we have developed 
a simple extension of the current LUBOEI algorithm that makes use of the 
production vertex distance between two mesons. 

Comparisons between the old model and the new one is made for vr-mesons 
according to exponential and Gaussian parameterizations of the BE correla- 
tion function. Qualitatively we note a relative suppression of the BE effects 
in multi-jet events, and especially for pairs of low momenta, with the new 
algorithm. This trend is in agreement with preliminary experimental data. 
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1 Introduction 



Like the Greeks thought of the world as consisting of indivisible parts called atoms, we 
today associate matter and forces with quantized objects simply called particles; or, in the 
mathematical language of Quantum Field Theory (QFT), particles are associated with 
excitations of fields. The theory that describes the subatomic world, as we know it today, 
is called the Standard Model. The simplicity of this model comes from the fact that there 
are only a handful of fundamental particles, such as quarks and leptons, which make up 
matter, and gauge bosons, which provide the forces or interactions. The complexity of our 
world comes from the fact that particles can take on different properties called quantum 
numbers, but more importantly they can be combined into larger units according to the 
rules set by the interactions. Combining three quarks one obtains baryons, combining a 
quark and an antiquark one obtains mesons. Mesons and baryons are members of a larger 
class of particles called hadrons. The combination of quarks into more complex particles 
cannot be made arbitrary, since they arc intimately controlled by the interactions; we will 
come back to this later. Of course, the everyday particles, that we are more familiar with, 
are not hadrons or leptons, but atoms that consists of protons, neutrons (baryons) and 
electrons (leptons), and molecules that are made up by several atoms. This last subject, 
however interesting it may be, is outside the scope of this thesis. 

Out of the four fundamental interactions only three are included in the Standard 
Model; the fourth, gravitation, still refuses to give in to all attempts of quantization 
that are based on standard gauge theories. However important in the everyday world, 
it is weaker than all the other forces in the microscopic world and can be neglected in 
the following. Electromagnetism and the strong force are treated in the theories QED 
(Quantum Electro-Dynamics) and QCD (Quantum Chromo-Dynamics) respectively. The 
weak force has no separate theory and was understood only when it was unified with 
electromagnetism, in what is known as the electroweak force. The theory of this unifica- 
tion is sometimes called QFD (Quantum Flavour Dynamics), or simply GWS (Glashow- 
Weinberg-Salam) . Strictly speaking one can regard the Standard Model as consisting of 
the theories QFD and QCD. 

The Standard Model is, as far as we know, an exact theory of phenomena accessible 
at current energies, but some of the mathematics contained within it is too complicated 
and cannot even be solved in a perturbative manner. In QED, the perturbative treatment 
has been so successful that physicists have claimed the electromagnetic force to be fully 
disclosed, whereas QCD is not so suitable for perturbative expansions, since it is a highly 
non-linear theory. For example, the running coupling constant Og, which is a dimensionless 
quantity that gives an effective strength of the theory, is larger than unity at low energies 
or, equivalently, at large distances. In the Feynman perturbative treatment, physical 
quantities are expanded in powers of ctg, which of course cannot be done for low-energy 
QCD. The origin of the perplexity of the strong force comes from the fact that the force 
carriers, called gluons, are massless and carry colour charge themselves, unlike the massive 
weak bosons and the (electrically) uncharged photons. These two properties mean that 
the gluons are self-interacting and can be produced without any significant cost of energy. 
In turn, this leads to a remarkable property of the strong force called confinement. Any 
quarks found in nature must exist in colour-charge neutral object (mesons and baryons); 
they are so to speak confined. Simply put, this property implies that there are no free 
quarks to be observed. To summarize, one can loosely say that QCD is a theory for 
quarks, but what we also need in order to understand the world is a theory for hadrons, 
the observable objects in our world. 
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Such theories can not be exact (nor unique) since they need to be mathematically sim- 
pler than QCD, hence they are called models. The more common ones are 'Lattice QCD', 
bag models and hadronization models. Here we will go into the domain of hadroniza- 
tion models. These are concerned with the process of creating hadrons out of their more 
fundamental building blocks, the quarks. Models of this kind are of great importance 
for experimental research, such as high-energy collisions at various accelerator facilities 
around the world. 

The model under consideration in this thesis is the Lund Hadronization Model, also 
called the Lund String Model or just the Lund Model jT]. The basic assumption of this 
model is that, because of the confinement property, a colour-field will mainly stretch out 
in one space dimension^, between a quark and an antiquark in the simplest case, forming a 
string-like object. The attractive force between the quarks will thereby become constant, 
denoted k, which is equivalent to a linear potential kt. We will come back to the more 
detailed concepts of this model later. 

Although the Lund Model is a semiclassical theory that only implements Quantum 
Mechanics in a probabilistic manner (not as a complex amplitude), it has been extremely 
successful and is widely used throughout the world by experimentalists as well as theorists. 
However, there are occasions where any simple probabilistic theory breaks down, i.e. when 
pure quantum mechanical effects manifest themselves. The Bose-Einstein effect is one of 
these phenomena that were never anticipated by classical physics. Quantum Mechanics 
requires that a class of particles called bosons (integer spin) need to have symmetrized 
wavefunctions under the exchange of identical particles, whereas the other part of the 
particle- kingdom, called /ermions (half-integer spin), has anti-symmetrical wavefunctions 
under the same type of exchange. A classical interpretation of this would be that of an 
attractive force for bosons, and a repulsive force for fermions, acting in both position and 
momentum space. 

Nevertheless, the Bose-Einstein effect has been implemented in the Lund Model in 
a pragmatic way by Lonnblad & Sjostrand 0. It is taken care of as a local shift in 
four-momentum between identical mesons according to a two-particle correlation function 
c((5), expressing the enhancement in probability to find a pair at a given Q, of the form: 

c{Q) = 1 + Ae-W^)" (1) 

Here A and R are phenomenological parameters, Q = a/— (pi — ^2)^ is the invariant 
four- momentum difference and t] is an integer that takes the values 1 or 2, to represent 
exponential or Gaussian shapes. This form is historically attributed to Hanbury Brown & 
Twiss [3^ who developed a way to measure star-disk radii using photon interferometry, but 
later the independent insight of Goldhaber et al. jl] led to its use in particle physics, as a 
way to explain the non-isotropic creation of pions in p-p annihilations. The Lund Model 
results, based on eq. (H)), has been verified by experiments for the case of 2-jets, but for 
higher jet-multiplicity it cannot fully reproduce what is seen. This is given by preliminary 
analyses of LEPl data which indicate that Bose-Einstein effects are overestimated 
by the LUBOEI algorithm in 4-jet events, when tuned to 2-jet events. More specifically 
this method seems to overestimate the Bose-Einstein effect for low energy hadrons. Also 
theoretical arguments gives that Bose-Einstein effect ought to depend on an 'actual' 
distance S rather than on some typical scale R. 

^This can be compared to the electric field that stretch out in three space dimensions, giving F ~ e^/r^ 
and V = — e^/r 
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In this thesis I shall try to resolve these problems by associating the constant R with 
the invariant four-position difference between hadron production vertices, called S (not 
to be confused with s, the squared centre of mass energy), up to a dimensionless constant 
k. For 2-jets there is an approximate correspondence between position and momentum 
which could explain why we can get by with one argument (Q) in the correlation function. 
For higher multiplicities of jets the complicated geometry of the string will break that 
correspondence. What is required is a two-parameter correlation function c(Q, S). We will 
show that such an approach will reduce the number of hadrons that have non-vanishing 
Bose-Einstein correlations in higher multiplicities of jets. The effects are not so big, 
however. Actually the results could also be seen as confirming that the simple ansatz (H)) 
is more meaningful than might have been expected. 

The outline of this thesis is as follows: Section 2 will give a review of the standard 
model, focusing on the properties of the strong force. In section 3 we introduce the 
features of the Lund Hadronization Model, of which the space-time picture and the Bose- 
Einstein model are our two main areas of interest. In section 4 we discuss an extension 
of the current Lund Bose-Einstein model and in section 5 the results of the extended 
model is presented. Finally in section 6 we summarize and give some ideas that might 
improve the simple extension for future work in this area. 



2 The Standard Model 

In this section I will describe the Standard Model in a simplified manner, by showing 
the known fundamental particles of the Universe and point out some interesting details. 
I will also try to explain the nature of the interactions by talking a little about groups. 
The following presentation is neither complete nor easy to follow for those having no prior 
knowledge, instead I advice [H] for a more complete treatment. 



2.1 Matter particles 

The matter particles, six leptons and six quarks, are distinguished by their flavours and 
are often arranged as doublets in three families or generations: 
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All the ordinary matter is made up by the first family. The families are similar in 
many ways but they differ in the particle masses. The third family has the highest masses 
and the first family has the lowest masses in each row. That is one reason why the families 
two and three are unstable, except for the neutrinos that cannot decay, only mix. 

The flavours of the matter particles could be said to be different values of one or several 
quantum numbers that represent mass eigenstates or weak force eigenstates, which are 
only approximately equal. The convention however, is to think of them as twelve different 
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particles, because the notion of mass is so important and so very different from other 
particle properties. 

Besides the flavour, quarks can have three colour charges: red (r), green (g) and blue 
(b). The leptons are colourless, or equally well white. Hadrons are also colourless in the 
sense that they have either all colours (r + g + b = white = colourless), or they have a 
colour and an anticolour. The concept of colour in particle physics is an extended analogy 
on Goethe's colour theory. 

Each particle has its corresponding antiparticle (the same symbol with a bar, or the 
same symbol with reversed sign for charged leptons), which is a particle with the same 
mass and total spin but with all other properties negated. Since the matter particles are 
fermions (spin-1/2) they have a helicity quantum number: right-handed particles (index 
R) have their spin parallel to the direction of motion and left-handed particles (index 
L) have their spin in the reverse direction. For massive particles, this last property is 
dependent on the observer, whereas particle/ antiparticle is a permanent property to our 
best knowledge (with some possible exceptions in the neutrino sector). 

2.2 Force particles 

The gauge bosons corresponding to the four fundamental forces are: 
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The graviton is not really predicted by the Standard Model but is included here for 
the sake of completeness. The photon is the carrier of the electromagnetic force and is 
its own antiparticle. The weak force has three gauge bosons W"*", W~ and Z'^, where the 
W's are each other's antiparticles and the Z° is its own antiparticle. Finally there are 
eight gluons carrying different combinations of colour charges, denoted by indices cc' in 
the table above. 

All the force particles are bosons and thus have integer spin; the graviton is a spin-2 
particle, whereas the others are spin-1 particles. 

The gauge bosons of the weak force are special since they have non-zero masses. In fact 
they are surprisingly heavy; the W's are around 86, and the Z° is about 97, times heavier 
than the Proton (938 MeV). Because of the Heisenberg uncertainty principle AiA£^ ~ h 
the range of the weak force is limited, thereby making it weak. 

What the charge of the weak force is might not be intuitively clear. In the table of 
the matter particles the parentheses are supposed to mark out weak iso-spinors, of which 
the single particles are merely projections. In analogy to the well-known theory of spin. 
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particles in the upper state of an iso-spinor has a third component of the weak charge 
being +1/2, and particles in the lower state has —1/2. This is however only true for 
left-handed particles; it is a strange fact that all right-handed particles have no iso-spin 
charge, i.e. they are singlets. 

Like the left-handed matter particles are doublets in an intrinsic iso-spin space, quarks 
arc triplets in a colour space. In this language different interactions corresponds to or- 
thogonal rotations in the intrinsic spaces; such rotations belong to mathematical groups 
called SU(n) if n > 1, and U(l) if the dimension equals 1. These are also called gauge 
groups, and in the Standard Model one usually give them an index that indicates which 
interaction the groups correspond to. 

As seen in the above table, the weak force, defined by the interaction particles 
W"*", W~, Z°), does not have an independent gauge group associated to it. The weak force 
is embedded in a larger gauge group called Uy(1) x SUl(2), which is the electroweak 
gauge group that also contains the electromagnetic group. 

2.3 Quantum Chromo-Dynamics 

In this thesis our main interest lies in the nature of hadronization, the process that creates 
hadrons out of the more fundamental quarks. This process is controlled by the strong 
force; it is included at the fundamental level in the QCD framework, but in reality it is 
less well understood. Nevertheless, I will state the basic equation of QCD and discuss 
some interesting properties of the strong force. 

2.3.1 The QCD Lagrangian 

The Lagrangian is in principle made out of the sum of Dirac^ Lagrangians of all known 
quarks and a Lagrangian term for the gluons. 

>Cqcd= Y1 qirT^^q-ltriG^'G,.) (2) 

g=u,d,c,s,t,b 

The last term is very similar to the Lagrangian of the electromagnetic field. The symbols 
for the quarks should be understood as wavcfunctions, or (spin-)vector fields, and the 
eight gluons are combined into one matrix field G^, all parameterized by space-time 
co-ordinates x'^ = {t, x, y, z). 

The covariant derivative, that substitutes the ordinary derivative, includes a gluon 
field term that causes the quarks and gluons to interact. 

1^u. = d^- igsG^ (3) 

The field strength tensor of the gluon field is a SU(3) generahzation of the electromagnetic 
field strength tensor; thus it is a 3 by 3 matrix in colour space. 

G^. = d^G, - d,G^ - i2g^[G,, G^] (4) 

The commutator term in G^u contains the physics that makes the strong force extraor- 
dinary compared to electromagnetism. It arises because SU(3) is a non-Abclian group 
(elements do not commute in general). It is an easy matter to show that, when the gluon 

■^The Dirac Lagrangian give rise to the Dirac equation, which is the correct relativistic generahzation 
of the Schrodinger wave equation for fermions 
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Figure 1: Different types of interactions in QCD: (a) 2 quarks and a gluon, (b) 3 gluons, 
(c) 4 gluons. 



field term of the Lagrangian is expanded, one obtains cubic and quartic terms in G^, that 
are not there for an Abelian group, such as electromagnetism. 

The field equations of the QCD Lagrangian are unsolvable. Especially complicated 
is the equation for the gluon field, which is cubic in (one order lower than the La- 
grangian). In practice, at high energies, one can use the Feynman rules to draw diagrams 
with vertices and propagators, and do perturbative calculations. Each term in the La- 
grangian with more than two fields is called an interaction term, which corresponds to a 
diagram where the particles meet in one vertex with a coupling strength as indicated by 
the pre-factors of the term in question. From the QCD Lagrangian it is clear that there 
exist three classes of interaction vertices. The first one is when a gluon couples to two 
quarks of the same kind. The 'same kind' can here also mean a quark and an antiquark 
of the same kind, depending on which direction you choose the time-axis in^. The second 
vertex is three gluons interacting, and the third one is four gluons interacting, as shown 
in Fig. m 

At low energies the Feynman rules are not applicable to QCD, instead one has to 
approach phenomenological rules to obtain a decent description of the physics. 

2.3.2 Confinement 

Just as charge neutral composite particles are energetically favoured by electromagnetism 
so are colourless particles favoured in the theory of QCD. Although it has not been 
proven from first principles, we have strong reasons to believe that the potential barrier 
to extract a free quark is infinitely high. The confinement property thus ensures that no 
quarks (colour triplets) or gluons (octets) can exist as free particles; they can only exist 
within mesons and baryons/antibaryons (all singlets). 

A meson consists of an antiquark with anticolour and a quark with colour. A baryon 
consists of three quarks, all with different colours. An antibaryon is then a combination 
of three differently (anti-)coloured antiquarks. The wave functions of the baryons/mesons 
has to be anti-symmetrized/symmetrized with respect to exchange of colours, so nothing 
precise can be said about what colour a single quark has. 

From an experimental point of view the confinement property is a hassle, because it 
hides the underlying physics from our eyes. But with modern high-energy accelerators 
some information on confined particles reveals itself in the form of jets. Jets are showers 
of collimated hadrons that arise after particle collisions. The formation of jets is a proof 

"^An antiparticle is equivalent to a particle travelling backwards in time. 
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that hadrons consist of quarks, and especially for n-jets with n larger than two it is the 
best proof of the existence of the gluon. 

For example, a collision that takes place between two protons primarily involves a few 
quarks and maybe some gluons. These particles will violently accelerate and change their 
direction of propagation, which is unmeasurable to us since this happens inside a confined 
'bubble'. But the bubble will break into pieces, and form hadrons that inherit the line of 
propagation from the accelerated quarks and gluons. Thus we see jets in our detectors. 

3 The Lund Hadronization Model 

Experimental as well as theoretical work has shown that the attractive force between two 
oppositely coloured particles is essentially independent of the distance at large separations. 
This fact is believed to have its origins in the non-linear gluon-gluon interactions in the 
surrounding colour field. The interactions cause the field lines to compactify into a flux 
tube (or string) of constant radius in order to minimize the action. One can associate a 
linear energy density k with the flux tube. More commonly we refer to k as the string 
tension, which is measured to have the pleasant value of ~ 1 GeV/fm. 

It is possible to use this fact to extract some physics, as it is done in the Lund 
Model. One could say that the basic ingredients to produce a feasible model that describes 
hadronization is: 

• A classical linear potential V{r) = nr. 

• A 'matter' string, on which gluons are distributed, and with a quark and an anti- 
quark attached at the two end points. 

• A semiclassical fragmentation model, approximating the correct Quantum Mechan- 
ical picture, such that it can be implemented with Monte Carlo methods. 

From these principles, in addition to several unmentioned, it is possible to consider 
some general classes of events that occur after particle collisions at high energies. In this 
study we restrict our attention to e+e" annihilation events. Then the simplest process 
is the creation of a quark-antiquark pair, e"'"e~ — > 7/Z*' qq. This is the two-parton 
case, which could for the moment also be called the two-jet case, but we wish to reserve 
this name for experimentally measured events; this difference in terminology will become 
clear later. The name parton historically applies to both quarks/antiquarks and gluons. 
Later it is worth discussing the 3-parton case, which will contribute with some interesting 
physics, that is used to extend our picture into a general n-parton case. 

3.1 The 2-parton case 

Assume we study an e~e"'" annihilating collision in its rest frame, with centre of mass 
energy E^ra) that creates a quark/antiquark-pair (qq). Assume further that the quarks we 
are discussing on the following pages are massless. The quarks will go out back-to-back 
(thus defining an x-axis) with the speed of hght (c = 1), spanning a gluon field string 
in-between. If we for a moment imagine that the string cannot break, the qq-pair will run 
out of energy when it reaches a separation of -Ecm/ and will be forced to turn back into 
a oscillating mode, as indicated by the dashed lines in the left part of Fig. |21 Since we 
have a linear potential and massless quarks, there is a direct way to convert between the 
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Figure 2: Left: A typical event that give rise to a qq-pair. Right: A zoomed in picture 
of the string as it fragments into hadrons. Note: space-time and energy-momentum 
co-ordinates are used interchangeably 

energy-momentum (in GeV) of a quark and the straight sections of its space-time world- 
line (in fm), with the help of k = 1 GeV/fm. Thus the 4-momentum and the 4-position 
can neatly be visualized in the same diagram. 

In reality, for the collisions studied in the Lund Model, Ecm ^ "^hadron for any 
metastable hadron. Thus the string will break up into smaller fractions even before the 
original quarks reach maximum separation. The breaking up is caused by the production 
of virtual qq-pairs that acquire 4-momentum from the string, allowing them to become 
real and screen off the colour of the two end-quarks. In principle a new qq-pairs could 
be produced anywhere in the string, but the breaks have to be correlated such that the 
hadron produced between two breaks has the correct physical mass and such that the 
total energy and momentum of the system is conserved in the fragmentation process. In 
the Lund Model this is taken care by an iterative Monte Carlo algorithm, that could be 
sketched in the following way. 

It is assumed that q^qj-pairs are created, where i, counted from right to left, runs from 
1 to n — 1, n being number of hadrons. Hadron i is then made up of the quarks q^qj-i, 
where q„ and qo are the two original quarks. 

If we let ko denote the 4-momenta of qo and fcg that of q„, then the first hadron, with 
4-momentum pi and mass mi, uses a fraction zi of ko when it is created. Zi is a random 
variable with probability density function 

f{z)cx-{l-zre~^ (5) 

where a and b are phenomenological constants; a ~ 1 and 6 0.5 GeV^^ . Also the hadron 
uses a fraction z[ out of fcg, where z[ is fixed by the restriction pf = {zik^ + z'^k'^Y = ml 
to: 

The second hadron (p2 and 1712) will have the remaining ki = {1 — Zi)kQ and k[ = 
{1 — z[)kQ at its disposal. These 4-momenta substitute ko and /cq above; other than that. 
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the reasoning is the same as for the first hadron. Thus one can define the remaining 
4-momentum ki, after i hadrons have been produced, to be: 

i 

h = n(i - ^^o^o (7) 

Here, the same equation apphes to k[ if one primes fco and Zj. 

The iterative procedure outhned above is easily continued any number of steps, so that 
every hadron (pj, mi) has its associated fractions Zi and z[. It can be continued as long as 
the remaining 4-momenta allows more hadrons to be produced; this gives us the number 
of hadrons n, which will vary from event to event. The constraints coming from overall 
energy and momentum conservation implies that not all Zi can be chosen independently 
according to eq. (0): the last two hadrons are necessarily fixed by the constraints. To 
minimize any non-physical effects of the choice of 'last two', fragmentation processes are 
allowed to begin in the two ends of the string and meet up at some point that is randomly 
selected. 

In the space-time picture (Fig. |2)) one can calculate the vertex where a q^qj-pair is 
created. If the vertices are counted from the upper right corner, then Wi = ki/n + (/cg — 
k'j)/K is the vector of the ith vertex. [Throughout this thesis we will consistently use 
the notation Wi for break-up vertices (qiqj-production points), and later t>j's will denote 
hadron production vertices.] From Fig. |21 it is clear that hadrons are represented by 
rectangles with areas equal to their mass-square, and with time-like diagonals equal to 
their 4-momenta Pi. When created, each hadron is connected to two neighbouring hadrons 
through the qq-vertices that define it. Of course, the hadronization picture used here is 
somewhat too ambitious, because the Heisenberg uncertainty principle does not allow us 
to speak simultaneously about the momentum and the position of the quarks inside a 
hadron. Thus the space-time picture should not be taken too literally; it is primarily 
used in the Lund Model as a way to deduce the 4-momenta of the hadrons in a particular 
event. Later we will try to expand the interpretation somewhat, such that one might 
speak of an approximate creation position of a hadron, which would be represented by 
some point inside the dashed rectangles in the right part of Fig. |21 

To make the model more realistic, a transverse momentum is added to each quark, 
with a Gaussian distribution p± ~ 300 MeV, in a manner such that it vanishes for each 
gg-pair originating from the same vertex. As a consequence we have to replace by 
the transverse mass = and all the 4-vectors used have to be projected onto 

the plane spanned by the two original 4-momenta k^ and k^, everywhere in the above 
equations and arguments. 

3.2 The 3-parton case 

If we have a three-particle system, e.g. created by e+e^ — — > qqg, it will in general 
spread out in two space dimensions (the xy-plane). In the Lund Model this is represented 
by a string with a kink on it. The kink/gluon carries 4-momentum, colours and spin, of 
which the first is central in the following discussion whereas the last property is neglected. 
The actual colours of the gluon are not very important, but for the sake of our discussion 
we may assign colours to the partons according toq = r, g = rb, q = b; this implies that 
the string flows from q to g and onto q. 

The kinematics of the string is more easily visualized if one boosts a three-parton sys- 
tem such that the momentum directions of the quark and antiquark is parallel respectively 
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antiparallel to the x-axis, and the gluon momentum is parallel to the y-axis. 

As seen in Fig. El the particles will travel in their own directions unaffected by what 
is going on at other point of the string. Meanwhile, they lose energy and momentum 
to the string being pulled out by them. Eventually the particles will run out of energy. 
Here, we assume that this happens to the gluon first, which is normally the case, but this 
does not do the argumentation less general. The gluon stops after travelling a distance 
corresponding to half its energy: since it has two string-ends attached to it, it loses energy 
twice as fast as a quark. As a consequence a third string region is formed that carries 
no momentum; it has only energy in this choice of reference frame. There are also two 
kinks formed when the gluon stops; they are however different from a gluon kink, since 
these two carry only information about the change of motion and do not have a local 
4-momentum concentration that is usually associated with particles. 

Later the quarks will also run out of energy, forcing them to stop. Then the virtual 
gluons in the string will finally catch up and make their presence felt. Since the string 
pieces closest to the quarks have a transverse momentum inherited from the gluon, the 
quarks start to move in the direction of the original gluon-momentum. The quarks will 
eat up a piece of the string and gain energy-momentum until they reach a point where the 
string is left behind again. The process repeats itself once more: a deceleration followed 
by a stop and change of path, followed by acceleration. Eventually all the three particles 
will meet in one single point translated away from the origin. This corresponds to a half 
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Figure 4: Projections of space-time sheet. The right picture is rotated 90 degrees in the 
x^-plane, with respect to the left picture. The arrows represent the parton world lines, 
the plain lines represent non-gluon kink world lines and the dashed line indicate a hidden 
part of the sheet. 

cycle of the gluon-quark-antiquark system at t = 7 in Fig. El The total momentum of 
the system exactly corresponds to the translated distance, which can easily be seen by 
adding up the original momenta that are indicated with dashed arrows (only half of the 
gluon momentum is showed). 

In this choice of reference frame it looks like a coincidence that the partons and the 
kinks travel in directions always parallel to the original momentum of the three particles, 
but this is of course true in all frames (for massless partons) . In fact the space-time world 
lines of the partons and the kinks can be built up by shorter lines that correspond to these 
three original 4-momentum vectors. The half cycle space-time sheet of the string will then 
look like Fig. |3] There are five different regions depending on which pair of four- vectors 
that spans them. They all have a diamond shape and can separately be boosted into their 
own rest frames, such that one obtains something that looks like the 2-parton case. 

As in the 2-parton case we consider only systems with high energy, which are the ones 
where the string model is valid. The 3-parton system will thus fragment into smaller 
pieces, corresponding to hadrons, at an early stage and will never evolve as far as Fig. |3] 
shows. The fragmentation process is essentially the same as in the 2-parton case for the 
two lower regions. In the middle region the available 4-momentum is not defined by the 
two vectors that span it, but by the remaining fractions of these vectors that is not used 
by hadrons in lower regions, thus reducing the number of hadrons that can produced 
with the fragmentation algorithm in this region. All fragmentation in later regions can 
be neglected in this treatment. 

When one becomes familiar with the geometry of Fig. 0] one can proceed by flattening 
out the first three regions, as shown in Fig. El so that an abstract but clear-cut picture 
emerges. This has to be done by cutting along the gluon 4-vector; hence it is mapped 
onto two vectors in this diagram. Also the 4- vectors are rescaled and put in perpendicular 
directions, corresponding to boosts into local rest frames. It is now easy to see the 
resemblance with the 2-parton case, and fragmentation is done as usual in each region, 
which is indicated in Fig. El The hadrons that cross over region boundaries are dealt with 
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Figure 5: An abstract picture of the 3-parton geometry. The area below the zigzag hne 
corresponds to the unfragmented string. The dashed hues indicate hadrons that are about 
to be formed. 




Figure 6: A 5-parton 'pyramid' that illustrates how an n-parton case might look in the 
abstract formalism. The four-vector pair that defines one region is shown as an example. 



by simply adding the two separate region contributions; together the regions contribute 
to the full mass-square and 4-momentum of this hadron. A main complication is that the 
z variable introduced in eq. © has to be generalized to allow for steps not only inside 
but also between regions. This can be done by relating z to the proper time of q^qj jZj, 
but we will not go into details here. 

It then becomes an easy task to generalize this discussion into an n-parton system, 
shown in Fig. IHlfor n = 5. Like the 3-parton case there will be more allowed regions than 
one might first expect. The number of regions equals the number of pairs of partons you 
can pick from the given n particles, disregarding the order in which they are picked. Thus 
there are nin — l)/2 available regions and they can be identified by the pair of parton 
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4- vectors that span them, i.e. {pi,Pj)- Fragmentation can occur in all the regions, but is 
less likely the further up in the 'pyramid' you are, and of course depends on whether the 
string has had time to fragment in lower regions, or did not have the time to do that. 

3.3 The LUBOEI algorithm 

There are several ways to introduce Bose-Einstein effects into a hadronization model. 
The most direct approaches are to define some simplified complex amplitudes, from wave 
functions or matrix-elements, that can be symmetrized for identical bosons before cal- 
culating the needed probabilities. According to Quantum Mechanics, these probabilities 
are proportional to the square of the absolute value of the matrix element, times the ap- 
propriate phase space factors. These models are often said to be of the global type since 
each event is associated with a weight that depends on the production of all the particles 
in the event. Global models has been tried out by people in Lund jH] and by others, but 
they all struggle with computational limitations and some inconsistencies that the models 
produce. 

The models of local type alters the correlation of 4-momenta between identical pairs 
of mesons within an event, and in contrast to the global types all generated events have 
unchanged weight = 1. The model worked out by Sjostrand & Lonnblad j2], sometimes 
called LUBOEI, is of this type. The local approach is more brute-force than the global 
one, but is computationally so much simpler that it has allowed more detailed studies. 

It has been found that a sensible way to account for Bose-Einstein effects is to intro- 
duce a two-particle^ correlation function c jH]. In its most simple form, which is useful 
for experimental purposes, it is a function of the four-momenta of the two particles. 

C(P1,P2) = p^^^,pf\ = 1 + Xf{pl,P2) (8) 

Here the numerator of the second expression is the double probability density to find 
particles with four- momenta pi and p2, given that the combined wave function is correctly 
symmetrized. The denominator is the product of two independent one-particle probability 
densities. In the limit where there is no Bose-Einstein effect this ratio ought to be unity. 
However, correlations from the jet structure of the event, energy-momentum conservation, 
etc., can in reality alter this significantly. In practice one therefore is more likely to 
define c(pi,p2) = Pbe^Pi, P2) / PnoBe{pi, P2) , where the numerator corresponds to the real 
world, with Bose-Einstein effects, and the denominator to a hypothesized one, identical 
in every respect except the need to perform a Bose-Einstein symmetization. A suitable 
parameterization is then the last expression in eq. (jH}, where A is an ad hoc parameter; 
naively A = 1 if / is normalized to f{pi = P2) = 1, but screening effects, such as when 
the two particles originates from different sources, will reduce the Bose-Einstein effects 
and give an effective A < 1. 

The function f{pi,P2) contains all the physics and should mainly depend on the invari- 
ant four- difference, Q = a/— (pi — ^2)^, of the momenta pi and p2. Goldhaber et al. jlj 
derived an approximate form, when considering a Gaussian shaped source, exp(— r^/2i?^), 
of radius R: 

fiQ) ^ e-«^«^ (9) 

^In the this section, whenever referring to 'particle', we will mean boson, or, more specifically, we will 
always talk about pairs of identical mesons. 
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In our studies we alternatively allow a correlation function where f{Q) is an exponential 
factor, instead of the Gaussian factor in eq. 0, and the final expression for c{Q) is the 
one given by eq. ((H). 

The experimental value of R is about 1 to 0.5 fm, or 5 to 2.5 GeV^^;thus the region 
where the Bose-Einstein effects are present is very small, Q < 0.2 — 0.4 GeV, compared 
to the dimensions of an extended string (~ 91 GeV at LEPl), and is comparable to the 
width of the string (or for Q, the transverse momenta of quarks). 

As a starting point we now consider uncorrelated particle production, uniformly in 
phase space, oc d^pi/Ei. A boost to the rest frame of a pair reveals that the relative 
momentum then is distributed according to 

d'Q QMQ 
— oc — = 10 

^pair v^4m2 + g2 

where m is the mass of a single meson. When Bose-Einstein effects are included, the 
expression in eq. ()10p has to be multiplied by the correlation function c{Q) to give the 
right distribution. A pragmatic way to include Bose-Einstein effects, when given a Q- 
distribution without such correlations, is to perform a local shift Q Q — S{Q) given by 
the equation 

r-Q q2^q rQ+S{Q) ^2^^ 



;i + Ae-(^«)^)— ^^^= / (11) 
y Am? + q'^ Jo y Am"^ + 

This may be solved approximately for small (5's, which yield 

Note that the algorithm does not change the number of particles, which implies that 
the effective c{Q) is less than unity for intermediate Q, contrary to the parameterization 
of eq. ()8I9|) . In the LUBOEI algorithm a shuffling procedure transforms the pairwise 
relative shift into a total shift for each hadron, owing to the net effect of all the pairs the 
hadron belongs to, such that the total momentum of the event is conserved. This also 
introduces some non-linear effects in the model, e.g. non-trivial three-particle correlations. 
Unfortunately, these types of operations cannot be done without violating either the 
conservation of energy or the conservation momentum of the total system. This is perhaps 
the biggest disadvantage with local models. The problem is dealt with by scaling all the 
hadron 3- vectors with a common factor in the rest frame of the event. However, there 
are other more sophisticated ways to shuffle the hadron 4-momentum, such as the one 
described in j2j, which could be argued to be more realistic, but for our purposes this 
simple procedure is enough. 



4 Identifying R with a space— time scalar S 

In the Goldhaber derivation of the correlation function, R is the radius of a spherical 
symmetric Gaussian source out of which pions emerge. At first glance the Lund String 
Model seems to lack such a source; instead the string is more of a cigar-shaped source, thus 
making it difficult to motivate why this correlation function should work well for 2-jets. It 
has been shown that the LUBOEI algorithm can reproduce experimental data with A = 
1.35 and R = 0.6 fm PU]; which tells us three things: First, the model is not fiawless, since 
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A should be at most unity in the interpretation given in section thus it is compensating 
for some physics that is still unaccounted for. Second, R is approximately the separation 
between neighbouring (neutral) and next to neighbouring (neutral & charged) pions at 
time of creation in the 2-parton picture. This is an interesting fact since pions are the 
hadrons with the highest multiplicities and thus give the dominating contribution to the 
Bose-Einstein effects. Third, R is comparable with width of the string (~ 0.7 fm) within 
errors; hence if R is to be interpreted as an effective source radius then this source is small 
enough to fit inside the apparent cigar-shaped source, which makes is possible to consider 
that it has a spherical symmetric distribution as a first approximation. The are however 
experimental data that indicate that the source has a certain elongation in the sting 
direction JT] as one might suspect, but we choose to ignore this small correction in this 
work. It should also be noted that physical considerations require that the Bose-Einstein 
correlation between well-separated hadrons is unmeasurably small, whichever source form 
is used; this explains the apparent misfit of the effective Bose-Einstein source and the 
cigar-shaped string region. 

If the effective source is not equal to the apparent source, defined by the string ge- 
ometry, then it is no longer necessary that the size of the effective source is constant 
for different pairs of hadrons within an event. The problem that Bose-Einstein effects 
for higher multiplicities of jets are not well described by the current LUBOEI algorithm 
indicate that the spatial structure of the string is important and that information about 
this should be invoked in the correlation function. The simplest way to achieve this is to 
parameterize c with two independent variables Q (kinetic) and S (spatial), where Q is the 
invariant 4-momentum difference as before and S, replacing R, could be something like the 
invariant 4-position difference of the hadron production vertices, S = a/— (^2 — ViY, see 
further discussion below. Keeping an expression very similar to the Q-parameterization 
of the correlation function, c{Q, S) becomes 

c(Q, S) = 1 + Ae-^*^'?^)" (13) 

where /c is a dimensionless constant, which is to be fitted to data. 

4.1 Production vertices of hadrons 

The experimental accuracy allows only a certain crude classification of production vertices 
(the points from which the particles seem to originate). Most particles appear to come 
from a common primary vertex, and only long-lived particles that have time to travel at 
least some millimetres give rise to resolved secondary vertices when they decay. Examples 
of the latter are D and B mesons, and A°, and particles like and vr^ that decay 
outside the detector. 

For theoretical purposes the experimental primary vertex is large enough to be resolved 
into separate productions regions for primary hadrons (hadrons that directly emerge out 
of the string fragmentation process) and secondary hadrons emerging from fast decays of 
the primary hadrons, e.g. p, K*, r] and higher excitations. The Heisenberg uncertainty 
principle prevents us from talking about a specific point where a hadron is produced with 
certainty, but it still allows us to speak about the mean value of a production point as if 
it was a classical quantity. Thus we can use the convention that a production vertex of 
a hadron is to be taken as the average of the area where the probability to produce the 
hadron in question is non-vanishing. 

Given that this area is approximately the same as the one that appears in the dashed 
region of Fig.[7| there are three natural points associated with a hadron production vertex. 
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Figure 7: The three candidates for a hadron production vertex. This picture represent 
the same 2-parton shown in Fig. |21 

An obvious choice, no. 1, is the point wliere the two constituent quarks meet for the first 
time. A second natural choice, no. 2, is the latest point which is casually connected to 
both surrounding qq-vertices. The third point, no. 3, is simply the average between the 
first and second point, which is the same as the average between the two qq-vertices, i.e. 
Vi = {wi + Wi^i)/2, where Wi for q^qj production and Vi for hadron i. By symmetry the 
production vertex should lie somewhere on the line that connects the first and the second 
point, and there is no reason to believe that a hadron is produced earlier than the second 
point or later than the first point. If the production vertex should be considered to be 
an average of the possible production points then the third point is in some favour. The 
quantum mechanical considerations in [H] also favour this choice. 

Fig. ISl is produced using P YTHIA [T2] , the Lund Model event generator, and shows 
that for low-Q pairs there is essentially no difference between S distributions of the three 
possible vertex definitions. It is mainly for pairs with low Q that the Bose-Einstein 
correlation is important according to eq. (fT!?|) . thus the physics will not crucially depend 
upon which choice we make. As a curiosity Fig. IHl shows the average of 5* as a function of 
Q in the fully allowed range. For Q larger than ~ i?cm/10 there is a clear divergence in 
the three vertex definitions, and at the extreme Q = E^m the different 5* are completely 
fixed with the values Ecm, 0, -Ecm/2 for points 1 to 3 respectively. 

From a theoretical point of view the third point seems more appealing, but the choice 
is not crucial for our discussion. As default we will chose the third point when talking 
about the associated production vertex for primary hadrons in the following pages. 

The fast decays of resonances, i.e. of p and K* in practice, might also contribute to 
Bose-Einstein effects, hence production vertices for the decay product must be deduced 
from simple principles. There are two ways of doing this: 

1. It can be argued that most of these decays are fast enough so that it is a good 
approximation to assume that the resonance decays immediately; in this sense the 
decay products are part of the primary hadrons and the resonance is merely an 
enhancement of the pair/triplet-mass distribution. The calculation of production 
vertices is then easily taken care of by inserting extra qq-vertices where appropri- 
ate inside the resonance production region, such that kinematical quantities are 
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Figure 8: Distribution of S = a/— (^^2 ^ ^i)^ for pairs of charged identical primary pions 
with Q < 0.5 GeV for a 2-parton event (only u-quarks, uii) at i?cm = 91 Gev. The three 
different cases correspond to the possible hadron production vertices shown in Fig[7| 



preserved. 

2. A more realistic picture would be obtained if the resonances where allowed to live 
for a proper time r, which is an exponentially distributed random variable with 
expectation value tq = l/F, tq being the lifetime of a particle and F its width. The 
decay vertex of a resonance, which in this simplified treatment could be taken as 
the production vertex of the daughter particles, is then given by: 

t'decay = h t^r (14) 

Here r stands for resonance and is the production vertex of the resonance. 



The production point assignment of a secondary hadron could be taken as an average of 
the possible decay states. For a decay resulting in two particles the two decay states (in 
1+1 dimensional space-time) are each others space reflections in the rest frame of the 
decay; thus the average production point must lie on the extrapolated world line of the 
decaying particle. The decay point, with method no. 2, is then a good approximation of 
the production point for both daughter particles. A similar argument hold in real space- 
time and for more than two decay products. This assignment would be unfortunate if any 
two of the decay products where identical mesons, since it would result in an overestimate 
of the Bose-Einstein correlation. As it turns out, none of the fast decaying resonances 
have allowed decays that would give this problem. 
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Figure 9: The mean (S) shows that there is an approximate correspondence between Q 
and S for charged identical pions. Here for the 2-parton case (-Ecm = 91 GeV, uu) and 
for three different definitions of the hadron production vertex. The error-bars give the 
one-sigma spread in each bin. 



4.2 A Lorentz invariant similar to the production vertex differ- 
ence 

When deahng with decays according to the second method (used as default) one runs into 
problems, since hadron vertices are no longer separated with strictly space-like distances. 
Quite often decay vertices end up inside the future light-cone of one or several of the 
primary hadrons. The ordinary invariant distance between production vertex (^ordinary = 
\/—{v2 — viY) cannot be the correct spatial scalar to use in eq. (fT!^ : firstly, S is badly 
defined for time-like distances, and secondly, the Gaussian form of the correlation function 
will blow up for large time-like distances. All this could be avoided if the first method of 
assigning production vertices for decay products is invoked instead of the second method. 
But that would be to avoid the problem by choosing a worse description of the physical 
process. Furthermore, in multiparton configurations, it can even happen that a primary 
hadron lies inside the future light cone of another primary hadron from a different string 
region. 

The conclusion drawn is that we have to define a new effective spatial invariant S^s that 
can be looked upon as the ordinary invariant 4-distance between two new effective vertices 
v[ and v'2, which in turn are functions of the production vertices as well as other relevant 
parameters. The two requirement that must be fulfilled is: Firstly, Ses and 5'ordinary 
must be approximately equal for pairs of primary hadrons originating from the same 
string region, at least when the exponent of the correlation function is small, kQS < 1. 
Secondly, Ses must necessarily be space-like, i.e. v[ and V2 must be simultaneous events 
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Figure 10: A space- like scalar 5'efr 
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in some reference frame. 

We have found it sensible to boost a pair of identical particles to their common rest 
frame and there define an effective Bose-Einstein distance 5'efr to be the invariant spatial 
distance between the particles at the time of creation of the last particle. This distance, 
illustrated in FigHUl then becomes 



where 1 is the early and 2 is the late particle. Here our choice of effective vertices were 
v'l = Vl + {{v2 - Vl) ■ {pi + P2))Pi/ ((Pi + P2) ■ Pi) (= ^^1 + (^2 - ti)pi/ El in pair rest frame) 
and V2 = V2- The new vertex v[ is then the point where particle 1 is found as particle 2 
is simultaneously created in their common rest frame. It is clear from the definition that 
S'cfr is always space-like. 

Fig. ^2 shows a scatter plot comparing the ordinary invariant distance 5'ordinary with 
the effective distance S'cfr for primary charged pion-pairs, for the 2-parton case and in the 
region of interest. The correspondence is satisfactory; a large portion of the points fall on 
the ideal diagonal, and the plot is symmetric for S'ordinary/efr < 1-5 fm, which is appealing 
because on average no bias will show up. At large distances the data shows (region not 
fully included in plot) that S'cfr > S'ordinary with increasing likelihood. 

From here on we shall only make use of the effective invariant distance defined in 
eq. (fT^ and hence we skip the index and refer to it simply as S. 

4.3 An algorithm to deduce primary production vertices 

Since the Pythia event generator does not keep in memory the primary production 
vertices, these have to be reconstructed from the output available. In this section we 
will show that the hadron 4-momentum, the rank of the hadrons and the original parton 
4-momenta will give us all the information required to calculate the needed production 
vertices. The parton 4-momenta will give us information about the string geometry; each 
parton pair spans a sub-string (or sub-sheet). In the Lund Model the notion of rank is 
used to specify the space/fiavour-order in which primary hadrons are produced. 

4.3.1 The 2-parton algorithm 

In the two-parton case the algorithm is simple: enumerate the primary hadrons from right 
to left with a rank i that runs from 1 to m. Call the qq- vertices Wi, where Wq = ko/n 
is the turning point of the original quark, and ko is the 4-momentum of the same, as in 




(15) 
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Figure 11: The approximate correspondence between ^ordinary and S^s for identical pri- 
mary charged pions in the 2-parton case (uii). 



section Then the first production vertex vi, as defined by point no. 3 in Fig. d is 
given by Vi = {wo+Wi)/2 = wo+Pi/2k, where pi is defined by the relation wi = Wq+Pi/k. 
This can be generalized to any hadron production vertex Vi with the recursive formulae: 

pi 

Wi = Wi^i H 

K 

wo = ^ (16) 

K 

Wi + Wi-i 
Vi = 

2 

i = 1 . . . m 

The translation vector pi between two neighbouring qq-vertices is the space-like diagonal 
of hadron i's production region, and is related to the hadron 4-momentum pi in the 
following way: 

~ Pi ' ko , , Pi ■ 

= T — U^o ~ 1 — ^7^0 (17) 

/Cg ■ Kg Kg ■ /Cg 

where fcg and /cq are the 4-momenta of the original quark and antiquark, respectively. We 
have in this equation assumed that /cg and fcg are light-like vectors, i.e. the quarks are 
massless. It can be shown, from conservation of the total system 4-momentum, that the 
recursion formulae give Wm = k'o/n, which is the correct answer. 

4.3.2 The n-parton algorithm 

For the ra-parton case the algorithm becomes much more technical, but the idea is the 
same as above. The physics discussion in subsequent sections can therefore be understood 
even if you choose to skip the rest of this section. 

The procedure of the 2-parton case must now be performed in separate regions spanned 
by two partons, but it cannot be done independently since hadrons may cross region 
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boundaries. Although there are n{n — l)/2 possible regions, the string fragmentation 
breakups can only occur in less than 2n — 3 of those regions, different for different events. 
This is due to the fact that fragmentation in an early region excludes fragmentation in 
some later regions, which are at time-like separations with respect to the first one. Hence 
a major complication, which one has to deal with in an n-parton algorithm, is to deduce 
in which regions the fragmentation (qq-production) actually occurred. To do this we need 
some definitions and conventions that are far from obvious. 

With the risk of confusing^ the reader we define the n-parton 'basis'-vectors to be 
ki, i = 1 . . .n, where ki is now the original quark 4-momentum and kn is the antiquark 
4-momentum, but where the intermediary k's are only half of the gluon momenta; and m 
is no longer the total number of hadrons but the rank of a single hadron. 

If we take all the hadrons to be mesons, then a hadron m is made up by a quark 
and an antiquark from two neighbouring breakups, which occur in regions {kj/,kii) and 
{kj//, kill) respectively. We then use the convention that the expression 'hadron m belongs 
to {kj, kiy means that the relations j' < j < j" and V <l < I" are fulfilled. An alternative 
interpretation of this expression is that the 4-momentum of hadron m has a component 
that can be written as a linear combination of kj and ki. With this convention the region 
with largest indices, j and Z, that a hadron can belong to will be the region where the 
production of the antiquark occurred, and the region with the lowest indices will be the 
one where the production of the quark occurred. 

Then we need to define a scalar indicator that gives the fraction of the basis-vector 
kj (or ki if j and / are interchanged) that is used by the m first hadrons. In order to 
make this quantity well defined one has to assume that we already know that hadron m 
belongs to {kj, ki). [From the onset we know this to hold for the region of the antiquark 
vertex of hadron m — 1.] Then z is given by: 

m j I ^ 

z{j,l,m) = (^Pi-^k,-^ki + ki + kj + k,^ 



=1 



kj • ki 



From the definition it might not be clear what the purpose of this indicator is, but it 
is intended to be used in the following way: if 



z{j, 1,171) > and z{l,j,m) < 1 =^ hadron m belongs to {kj,ki) 

z{j,l,m) < and z{l,j,m) < 1 =^ hadron m belongs to {kj^i,ki) 

z{j,l,m) > and z{l,j,m) > 1 =^ hadron m belongs to {kj,ki^i) 

z{j,l,m) < and z{l,j,m) > 1 =^ hadron m belongs to {kj^i,ki+i] 



(19) 



In either of the three last cases, the procedure is to be iterated from the new region defined 
on the right hand side. If the indicator gives the first case one knows that the antiquark 
of hadron m was produced in the region defined by this case; thus implying that hadron 
m + 1 belongs to this region, and the procedure can be continued with this hadron. We 
have here assumed that 1 < j < / < (in the definition of z this is not assumed) since this 
will give us the correct number of regions. This implies that some of the regions on the 
right hand side will not exist, but with the physical picture described so far the left hand 
side criteria will never be fulfilled unless the right hand side statements make sense. Thus 
with the 2;-indicator it is possible for every hadron to iteratively deduce which regions it 
belongs to, which in turn will give which region a specific qq-vertex {Wm) belongs to. 



®The notation in this section should be looked upon as being independent of the one used in section lTTI 
although some parts are reused 
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Figure 12: The r^-distribution (where r is the proper time of hadron production vertex) 
for the ideal algorithm and the r-maximizing algorithm in a 3-parton 'Mercedes' event 
with = 91 GeV (uu). The distribution 5*^ here refers to the usual invariant distance 
between the two production vertices calculated by the two algorithms for every primary 
hadron. Note 1: a delta function contribution at the origin, corresponding to identical 
vertices, is not shown in the S"^ curve. Note 2: The apparent negative r^'s in the r-max 
algorithm is an effect caused by the drawing routine and the choice of bin-size, not a real 
result. 



With the definition of rank the first hadron will always belong to the first region 
(right-most region in Fig. Ej) labelled by {ki,k2)- The vertices are then calculated with 
formulae that can be made recursive like to those of eq. ()16p . but are better displayed in 
a direct form: 



z{j,l,m)^ + y] — 

K K ^ — ^ K 

i=j 



Wo = - (20) 

K 



2 

m = 1 . . . tl primary hadrons 

Here j < I are chosen so that Wm, which is the antiquark production vertex of hadron m, 
belongs to region [kj, ki); thus j and I are functions of m, determined with the aid of the 
indicator z. As before Vm is the m'th hadron production vertex. 

With an ideal n-parton system, such as described in section the vertex algorithm 
should be perfectly flawless. In practice, when using an event generator many further 
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Figure 13: As Fig. [T21 but with a physical mixture of n-parton events at Ecm = 91 GeV 
(using PYEEVT, only uu). 



physical aspects are included, that might or might not be important for our study, but 
cannot simply be turned off and will inevitably give rise to complications. It was discov- 
ered when using Pythia that for n > 2 our algorithm, with the strict criteria in eq. (fT^ . 
could not always choose the correct region to place a hadron in. A small fraction of 
the hadrons were assigned production vertices with a negative proper time, which is not 
acceptable. This can be related to ambiguities in the hadronization algorithm itself [Zj, 
where e.g. transverse momentum fluctuations can lead to contradictory assignments of 
production region. A pragmatic solution to this problem was to evaluate for each hadron 
several possible production vertices, corresponding to every possible region spanned by 
two partons, and then to pick the vertex with largest proper time. In principle you evalu- 
ate the Wm and Vm for every region {kj, ki) with the formulae in eq. (j2U|) . and then compare 
for each m all the = given by different j and /. A further minor approximation, to 
reduce the computing time as well as the code, resulted in that only the two lowest levels 
of regions (1 < / — j < 2) were evaluated in the process. The result was not a bias for high 
proper times as one might expect; instead the modified algorithm and the ideal algorithm 
produced in the 3-parton case a very similar distribution of the proper time, shown in 
Fig. ^1 with the exception that the former was liberated from a tail of negative values. 
For a physical mixture of n-parton events, there is a small but clear discrepancy between 
the definitions for the majority of data points, shown in Fig. ^| This is caused by the 
large number of soft gluons, which make the space-time picture fuzzy, thus undermining 
the validity of the ideal algorithm. Nevertheless, the two algorithms can be considered as 
being equally good (or bad) in determining the correct vertices, and the discrepancy can 
be interpreted as a measure of the uncertainty in the physical understanding. Henceforth, 
we choose to use only the modified algorithm and avoid any negative proper times. Note, 
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Figure 14: The phase-space of vr^'-pairs in 2-parton events (uu) at = 91 GeV. Right 
frame is a blow-up of the interesting region. 

4 




Figure 15: As in Fig. El but for vr^-pairs and vr -pairs. The number of events plotted is 
the same as used in Fig. HH 



even though we have assumed all quarks to be massless in the recipe of this algorithm as 
presented here, in the algorithm actually used to produce the results we also accounted 
for heavy quarks by finding the light-like vectors that span a given parton region, with 
r = at the primary production vertex. 

4.4 The shift (5 as a two- parameter function 

In order to invoke the Bose-Einstein correlation given by eq. ()1H|) we need to know how 
the phase-space looks like in the parameters Q and S. It is interesting to see the phase- 
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space of the basic event, the 2-parton event, even though we will not make further use 
of it. This is shown by the PvTHlA-generated scatter plots in Fig. El and Fig. for 
the two most common mesons. At first the Q-S dependence is not as clear as one might 
expect, but it is obvious that the phase space is restricted into a triangular area, of a size 
that depends on the system rest energy Ecm, and a shape that depends on the definition 
of the hadron production vertex. It is mainly the phase-space restrictions and the fact 
that low Q and 5* seems favoured that give rise to the behaviour of the mean (S), as a 
function of Q, shown in Fig. El The phase-spaces of 7r°-pairs and vr^-pairs are similar at 
large scales, but there is one difference that might become very important for their Bose- 
Einstein behaviour. Neutral pions pairs are closely clustered along the lower bound in S, 
which corresponds to nearest neighbour pairs, whereas charged pions are forbidden to be 
produced in such pairs. The same effect shows up in the less interesting upper bound of 
S: the flavour-order in which hadrons are produced in the 2-parton case prohibits charged 
identical particles from being produced at the two string ends. 

For multiparton events the phase-space will not look as simple as above, instead it 
will more and more take on a familiar shape coming from randomly distributed pairs in 
real and momentum space. Thus we have to continue our discussion assuming a small 
correlation in Q and 5*. 

If we wish to follow the path taken in section 13.31 in a more general manner we would 
look for a shift {6q,6s) in both variables. This would correspond to a situation where 
we were given a distribution in {Q,S), in a world without Bose-Einstein effects, and 
then were asked how this distribution with a minimal shift^ {Q,S) {Q,S) — {6q,6s) 
would be turned into the correct Bose-Einstein distribution. In the situation where the 
difference in momentum and position for identical pairs of bosons are isotropically and 
homogeneously distributed throughout both spaces the ratio between the correct and the 
non-Bose-Einstein distribution is given by the correlation function in eq. ()13p . according 
to our assumption. 

A two-parameter shift would only work if Q and S are independent variables, however, 
which is not entirely true. A shift in Q would induce an indirect shift in 5* and vice versa, 
hence this two-parameter shift is problematic. We have therefore chosen only to alter the 
Q-distribution, at the expense of not performing a minimal shift, i.e. one might choose 
a larger 6q instead of smaller {6q,6s)- Of course there is still an indirect shift in S, that 
we will ignore. That is, given an original event without Bose-Einstein effects, the S, 
Q and 6q are evaluated, but the procedure is not iterated to provide updated S values. 
Furthermore, Q is a measurable quantity, which is not true for 5*. 

We do not know for sure how the shift of momenta from Bose-Einstein effects are 
correlated with shifts in the space-time picture, so, in order to proceed, we make the crude 
assumption that S is constant under a shift in Q, thus giving a constant contribution to 
the phase-space density, which factors out. Then nothing changes from the derivation of 
S{Q) in section ISSl except for the parameterization of the correlation function, and one 
obtains a shift in the Q-variable given by: 



The integral is not easily solved in an analytically useful form. Instead we had better 
perform a numerical integration in a computer efficient way whenever implementing it. 




(21) 



minimal shift is a shift that alters the physical picture (string geometry etc.) the least, but still 
gives the correct distribution 
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5 Results 



We have evaluated the algorithms of the proposed new Bose-Einstein model using the 
event generator Pythia. Since the particular choice of phenomenological Bose-Einstein 
model ultimately depends on experimental comparisons our objective is to simulate a 
reahstic experimental situation and to produce data that are as close as possible to the 
one obtained in the real world. But since relevant experimental data are scarce, and since 
comparisons are best performed by experimentalists, we will only qualitatively compare 
the new [Q, S')-dependent algorithm with the current Lund model Q-dependent Bose- 
Einstein algorithm. To simplify the discussion we choose to denote the two models with 
BEqs and BEq, respectively, and NoBE is used to denote a reference null- model, which 
is the Lund Model without Bose-Einstein effects. 

The BEq-version used is the one with global energy compensation, which is referred 
to as BEq in jT2]. In the BEqs algorithm we use the same kind of energy compensation. 
Whenever comparing BEq and BEqs we will always set the switches equivalently, so that 
the models only differ by the S'-parameterization in the correlation function, e.g. when 
comparing exponential models we setrj = 1 and MSTJ(51) = 1, when comparing Gaussian 
models we set rj = 2 and MSTJ(51) = 2. Note: When developing the BEqs naodel we 
had the idea that S oc Q for small Q's, thus one might want to compare BEQs(exp) with 
BEQ(Gauss), since the the two correlation functions show similar Q-dependence. But, 
in fact, S approaches asymptotically a constant (~ 1 GeV) in the low-Q limit, and the 
above stated comparisons are more relevant. 

We will only discuss the lightest and most abundant meson, the pion, of which like- 
charged pion-pairs (tt+tt^ + tt^tt^) have the same Bose-Einstein behaviour, which is 
slightly different from the one shown by neutral pion-pairs (vr^vr^). Since charged pi- 
ous are more often used in experiments, being easier to detect and measure due to their 
charge, we will focus on the behaviour of these particles, and the data from neutral pions 
will be presented when the results are interesting enough. 

In the analysis we used the subroutine PYEEVT that generates random ra'-parton 
events according to realistic e"'"e~ annihilations at LEPl, at Ecm = 91 GeV. These were 
classified to n-jets using the subroutine PYCLUS (Scaled Durham distance: MSTU(46) 
= 6, Maximum joining distance: y = PARU(45) = 0.005 [12j) such that the variable 
number n was obtained. For statistical reasons all events with n > 6 were added to the 
5-jet events. 

Here it becomes obvious why we choose to make the distinction between an n'-parton 
event and an n-jet. The number of partons n' in an event is not an experimentally 
measurable quantity; furthermore, in principle, but not manifestly included in PYEEVT, 
there are an infinite number of partons in the limit where the gluon energy or emission 
angle goes to zero. For practical purposes all one can do is to define a suitable experimental 
criterion, as done above, that allows one to deduce the number of distinguishable jets n 
out of a larger number of indistinguishable partons n'. 

For these simulated LEPl-events the Q-distribution dN/dQ, N being number of pairs 
in an event, was collected for charged pion-pairs and neutral pion-pairs. This also includes 
pairs where one or both come from decays of long-lived particles and thus cannot display 
Bose-Einstein effects. The familiar Bose-Einstein ratio BE/NoBE, in the Q-distribution, 
was then obtained by: 

dA^BE ^ dNBE/dQ 

diVNoBE diVNoBE/dQ 
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Figure 16: The BEqs/NoBE and BEq/NoBE ratios in the Q-distribution for charged 
pion-pairs. Exponential models are used. 
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Figure 17: As in Fig. but for neutral pion-pairs. Exponential models are used. 
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Figure 18: The BEqs/BEq ratios in the Q-distribution for charged pion-pairs. As guid- 
ance for the eye the constant 1 is plotted as well. Exponential models are used. 



This ratio, which should resemble the correlation function c{Q) , is shown in Fig. and 
Fig- El for the two exponential models and for different jet multiplicities. Note that the 
number of n-jet events differ between the BE and NoBE alternatives, and that therefore 
the dN/ dQ distributions are normalized per event before taking ratios. 

Considering charged pions and exponential models. Fig. ^1 shows that the Bose- 
Einstein effect is stronger for higher jet multiplicities in both models, but that the new 
BEqs model suppresses the effects relatively to the BEq model. Since the BEq model, 
as far as known, reproduce experimental data in the 2-jet case we have normalized the 
BEqs model to give the same behaviour in this case for charged pions, by setting the free 
parameter k in the correlation function (jl3|) to = 1/11 (for k = 1 GeV/fm and recalling 
that h = 0.197327 GeV-fm). The value of k is relatively small, one would perhaps expect 
a value close to unity, but it might reflect the fact that in our model we have only used 
the average position of a hadron at time of creation, not the hadron wavefunction. The 
wavefunction is of course spread out in space and the effects of a symmetrization is most 
notable at points in-between particles, thus corresponding to a smaller S than we have 
used. The value of k is also dependent on other parameters, in particular R since we have 
fitted the BEqs model to the BEq model, which depends on R. The relation is that R 
and k should be approximately proportional. The relevant parameters used throughout 
the analysis are: the constant R in BEq correlation function (P) is set to i? = 0.6 fm; A 
in both correlation functions (jlll3|) is set to A = 1.35. 

For the neutral pions the Bose-Einstein correlation is even more suppressed than 
for charged pions in both models. This somewhat unexpected behaviour has a simple 
explanation in the decays of resonances that produce a multitude of uncorrelated neutral 
pions that will dilute the real effect in our way of presenting data. 
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The contrast is small between the two models in Fig. ^1 and Fig. El thus, instead we 
make use of a double ratio in order to enhance any differences: 

dA^BEqs ^ dA^BEgs/dA^NoBE ^^3) 
dA^BEq diYBEq/dA^NoBE 

This is shown in Fig. where one for charged pions can see a relative suppression in the 
Bose-Einstein correlation of the BEqs model. For 4-jets there is a 2-3 percent suppression 
and for 5-jets and higher there is a 4-5 percent suppression depending on the choice of 
Q-region. In the case of neutral pions, Fig. 1191 show a similar behaviour between different 
jets, but there is also an overall shift upwards. Note that the general level of Bose-Einstein 
effects always can be adjusted by a parameter retuning, so it is the trend with increasing 
number of jets that is the interesting feature. 

We have so far discussed the comparison between the two models with exponential 
correlation functions, but, in fact, quantitatively the results obtained with Gaussian cor- 
relation functions are very similar. In Fig. 1201 the BE/NoBE ratios and in Fig. 1^ the 
BEqs /BEq ratios are given for charged pions. The shape of the curves in this last plot is 
somewhat different from the equivalent Fig. for exponential models. But the 2-jet fit 
is consistent with an approximative value of = 1/9 (used in all Gaussian plots), which 
is slightly larger than the one obtained from the exponential fit. Thus it seems that the 
Gaussian parameterization in itself gives stronger Bose-Einstein correlation. Which of 
the two different types of shapes that is more true is ultimately an experimental question. 

For neutral pions with Gaussian models the result is, as expected, equivalent to earlier 
result: the Be/NoBE ratios are very similar to Fig. [T7I fneutral pions, exponential), and 
the BEqs/BEq ratios are very similar to Fig. 1711 (charged pions, Gaussian) with an overall 
increase (above unity), which corresponds to stronger Bose-Einstein correlation with the 
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Figure 20: As in Fig. but for Gaussian models and charged pions. 
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Figure 21: As in Fig. ^Ibut for Gaussian models and charged pions. 
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Figure 22: The BEqs/BEq ratios in the Q-distribution for charged pion-pairs in the 2- 
jet case (uu). The intervals give the relevant i?pair-bins, whereas the curve labelled by 
BEqs/BEq gives the total ratio for 2-jets, which also occurs in Fig. UHl but there for all 
quark flavours. Exponential models are used. 



new model. This we have already observed for exponential models in the neutral and 
charged pions plots. 

Since it is sometimes possible in experiment to distinguish e"'"e~ — > Z*^ — uii/dd/ss 
events from heavy quark events, e"'"e~ ^ Z° — > cc/bb, we have also looked at the for- 
mer events specifically. No extra plots will be presented since the result does not differ 
substantially, but one can generally say that the Bose-Einstein effect in this case is less 
diluted. There is approximately a 15 percent increase in all BE/NoBE ratios, but for the 
BEqs/BEq ratios the change is not significant. 

It should be noted that for all plots with the ratio BEqs/BEq, i.e. Fig. UHl Fig. [1^1 
and Fig. 1211 the ratios do not approach unity outside the range included in plots (Q > 1). 
Most of the ratios take some constant value in the interval [1.00,1.02] for Q > 2. However, 
it is not required that this ratio should become unity in an event generator, since for both 
models dN/dQ approach zero very fast in this range, and small effects are blown up. 

5.1 The effect of the BEqs algorithm in different energy bins 

The smallness of the BE-suppression for charged pions obtained with BEqs is impressive 
but discomforting: impressive that so stable, i.e. the naive BEq is not that bad, but 
discomforting since it becomes hard to tell if the suppression is significant enough. This 
leads us into further investigations of the actions of this new algorithm. A critical test 
is whether or not the new algorithm fulfils the following intuitive insight. For higher 
multiplicities, say 4-jets, the string regions are relatively well separated in space, but the 
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Figure 23: As in Fig. but for 3-jets. 
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Figure 24: As in Fig. |22lbut for 4-jets. 
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Figure 25: As in Fig. |221but for 5-jets and higher orders. Note: Fluctuations are statistical, 
not caused by model. 



central parts of each region holds relatively little of the parton momentum, thus these 
parts are closer in momentum space. Hadrons coming from these central parts have thus 
low Q and high which should enhance differences in a comparison between BEq and 
BEq, if the new algorithm satisfy our expectations. More specifically, we want to compare 
meson-pairs that inherit a small portion of the parton momenta, i.e. mesons in the centre 
of a 2-parton region, to pairs that inherit a large portion of the parton momenta, i.e. 
mesons closer to the partons in the space-time picture. 

In order to distinguish such pairs we choose criteria that are also useful for experi- 
mentalists. For two mesons with energy Ei and E2 in the total system rest frame, we 
define -Epair = -Ei + -^2- Then we create four energy-bins by cutting the allowed £^pair 
range into appropriate regions, given in Fig. |221 to Fig. |2H1 These plots show the double 
ratio BEqs/BEq for each separate i?pair-bin. 

By splitting the data into four regions like this we will introduce some unavoidable 
correlations, since Q and -Epair are not independent quantities. But the effect will be small 
unless Q is near a cut in the i?pair-range. It should be noted that the relative number of 
pairs in an ii^pair-bin vary from 2-jet to 5-jet, thus influencing the relative behaviour of 
the plots when comparing different jets in Fig. |221to Fig. |231 

If pairs with low-£'pair are interpreted as hadrons that have been created in the central 
parts of parton-regions, then it is clear that BEqs model will strongly suppress Bose- 
Einstein correlations from these pairs, e.g. 5 percent for 2-jets and 10 percent for 5-jets 
and higher. But, oppositely, high-energetic pairs, corresponding to hadrons created near 
region boundaries, will with the same model contribute with an intensified Bose-Einstein 
effect of similar strength. However, the net effect will always be a suppression for n-jets 
with n > 2, as already noted. Since the BEqs algorithm is roughly tuned to have the 
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same overall dN/dQ distribution as BEq, the supression of Bose-Einstein effects at small 
impair in the former then leads to an enhancement at large -Epair- Had the tune been for 
agreement at large -Epair, the drop at small i^^pair would have been even more dramatic, 
up to 20 percent for 4-jet events. 

We present only plots for charged pions and exponential models, with initial u-quarks, 
i.e. only e"'"e~ uu events (this increases i?pair-bin spread by approximately 

60 percent compared to the full mixture of quark flavours). Neutral pions-pairs behave 
essentially the same way as charged do in Fig. |221to Fig. 123 Of course, there is a general 
shift upwards in the curves corresponding to these particles, which can be seen already 
in Fig. ^1 versus Fig. ^1 In addition, for neutral pions there is a smaller spread between 
the curves coming from different i?pair-bins. The spread is about 85 percent of the spread 
seen for charged pions. However, this last effect is probably caused by the dilution effect 
from decays, mentioned earlier, that is stronger for neutral pions compared to charged 
pions. 

Also when comparing the Gaussian models the result is the same as obtained with 
exponential models. It is always the low- impair bin that drops down by a few percent 
compared to collective curve, and high- impair rise by approximately the same percentage. 
The spread in ii^pair-bins is approximately 75 percent and 60 percent for charged and 
neutral pions, respectively, compared to the spread seen for charged pion with exponential 
models. 

6 Summary & Outlook 

There has been an indication that the LUBOEI algorithm overestimates the Bose-Einstein 
effects for 4-jet events, when the parameters R and A is tuned for 2-jets |l5|. We have reason 
to believe that this phenomenon has to do with the fact that the space-time picture is 
neglected in the correlation function c{Q). This led us to try out the most simple extension 
of this model, by incorporating the space-time picture in one variable S, corresponding 
to an invariant separation of production vertices, thereby giving a correlation function 
ciQ,S). 

The Lund Model space-time picture used for string fragmentation provides natural 
hadron production regions, where a single hadron may be associated with a point that 
represents the hadron production vertex. 

Including decays of primarily produced resonances and multi-parton events showed 
that the quantity S cannot simply be the ordinary 4-vector distance of two hadron pro- 
duction vertices, since there can always be vertices that are at time-like separations as 
well as space-like. Instead a possible choice is that S is the 3-space distance between the 
two hadrons in the pair rest frame when the late particle is produced. 

The structure of the Lund Model string allows an algorithm to untangle all the hadron 
production vertices from information given by hadron rank, hadron 4-momenta and parton 
4-momenta. Unfortunately ambiguities in the hadronization algorithm occasionally gives 
hadron 4-momenta that cannot consistently be used for assigning production vertices with 
this algorithm. An easy but surprisingly functional solution to this problem is to evaluate 
a hadron production vertex in all possible parton regions and then choose the vertex that 
maximizes its proper time. 

The effect of the new correlation function depends on the phase-space given in vari- 
ables Q and 5", thus the result will depend on how we choose to perform local shifts in 
these variables. Due to uncertainties we made the crude assumption that the shift is in- 
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dependent of the 5" phase-space density and that the direct shift is only to be performed 
in the observable Q. 

The new BEqs algorithm then proved to possess many desirable properties. For 
charged pions the Bose-Einstein effect was suppressed in 4-jets and 5-jets compared to 
the BEq algorithm, whereas 2-jets and 3-jets gave comparable results. This is qualita- 
tively what is found in the preliminary LEPl-results of [Sj. For neutral pions the result 
was also an enhancement of the Bose-Einstein effect in general. However, in order for 
them to produce the same result for 2-jets the parameter R has to be tuned separately 
for charged and neutral pions, with the relation i?^± > R^^o . And indeed this is consistent 
with experimental results [TT]. 

In the new model we introduced a free parameter fc, with the rough value 1/11, that 
was fixed by fitting the BEqs 2-jet curve for charged pions to the equivalent BEq curve. 

We also found that the new algorithm acts differently depending on the energy of the 
pion-pair, -Epair- The Bose-Einstein effect is relatively suppressed in the low- impair region, 
which accounts for the majority of pairs, but enhanced in the high-ii^pair regions compared 
to the results given by BEq. 

The results so far are very interesting and we suggest that further work could be done 
by considering: 

• Study the effects of the new algorithm in e"'"e~ W^W~ events, which give two 
quark-ant iquark systems. These events have lately been given much attention in 
experimental studies, and they are very important indicators that show if the BEqs 
algorithm is on the right track. The issue is whether or not there is an inter-string 
Bose-Einstein effect. If there is not, then the two string pieces are independent 
systems (no cross-talk), if there is, then it cannot be modelled by a simple Q- 
dependent model, the spatial inter-string distance must show up somehow in a 
realistic parameterization. 

• Different parameterization of the correlation function, such as using an 4-product 
between a spatial quantity and a momentum quantity in the exponent in eq. (fT^ . 

• Modelling an elongation of the source in the same correlation function as suggested 
by experimental data [TT]. We have not studied it, but the BEqs model might 
already have such properties. 

• Study the phase-space density dependence on 5* and take local shifts in both Q and 5* 
into considerations when producing a Bose-Einstein influenced (Q, S')-distributions. 

• Study the role of different ways to include secondary decays, i.e. method 1 versus 2 
in section 14.11 and other possibilities. 
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